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In this paper we investigate the time evolution of the energy density for a real massless scalar 
field in a two-dimensional spacetime, inside a non-static cavity, taking as basis the exact numerical 
approach purposed by Cole and Schieve. Considering Neumann and Dirichlet boundary conditions, 
we investigate the following initial states of the field: vacuum, thermal state, the coherent state and 
the Schrodinger cat state. 
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I. INTRODUCTION 

Moore [l[ , in the context of a real massless scalar field inside a non-static cavity, obtained an exact formula for 
the expected value of the energy-momentum tensor, given in terms of a functional equation, usually called Moore's 
equation. For this equation there is no general technique of analytical solution, but exact analytical solutions for 
particular movements of the boundary an d also approximate analytical solutions (see, for instance 

have been obtained. On the other hand, Cole and Schieve [7f proposed a numerical method to solve exactly the Moore 
equation for a general law of motion of the boundary. Moreover, the dynamical Casimir effect has been investigated 
for a thermal bath as the initial field state in a non-static cavity 0, d, 0, [l(| El, and also for a coherent state 
[ToL [l2l [l3L HH . Squeezed states have also been considered • 

In the present paper, taking in account a technique of calculation based on the Cole-Schieve approach 0, fl5| . 
we examine the energy density in a non-static cavity for the following initial states of the field: vacuum, thermal 
state, coherent state and the Schrodinger cat state. The Cole-Schieve approach is appropriate, since we consider a 
certain non-oscillatory law of motion, with large displacements and relativistic velocities. In this sense, perturbative 
techniques requesting non-relativistic oscillating movements with small amplitudes are not suitable in this case. The 
article is organized as follows: in Sec. 2 we show the general field solution. In Sec. 3 we investigate the non-static 
cavity for a vacuum and a thermal initial field states. In Sec. 4 we analyze the coherent and Schrodinger cat states 
as the initial field states. In Sec. 5 we make some final comments. 



II. GENERAL FIELD SOLUTION 
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Let us start considering the field satisfying the Klein-Gordon (we assume throughout this paper h — c — ks = 1) 
and obeying conditions imposed at the static boundary located at x — 0, and also at the moving boundary's position 
at x = L(t), where x — L(t) is a prescribed law for the moving boundary and L(t < 0) = Lq, with Lq being 
the length of the cavity in the static situation. We consider four types of boundary conditions. The Dirichlct- 
Neumann (DN) boundary condition imposes Dirichlet condition at the static boundary, whereas the space derivative 
of the field taken in the instantaneously co-moving Lorentz frame vanishing (Neumann condition) at the moving 
boundary's position: We also consider: Dirichlet-Dirichlet (DD), Neumann-Neumann (NN) and Neumann-Dirichlet 
(ND) boundary conditions. A general solution of the wave equation can be written as: 



where the field modes ip n (t,x) are given by 



E 

=1-2/3 



a n ip n (t, x) + B.c. 
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with t/>(°) = [R(v) + R{u)] jl [Hj], u = t — x, v = t + x, and R satisfying the functional equation R[t + L(t)\ 
R[t — Lit)] = 2 , which is the Moore equation. The operators A and B satisfy the commutation rules A, B 



1, 

0. The NN solution is recovered for A = 7 = 1 and (3 = 0. The other three cases are recovered if 
A = and: (3 — and 7 = i for the DD case; (3 = 1/2 and 7 = 1 for the DN case; (3=1/2 and 7 = 1 for the ND case. 
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FIG. 1: (a) Moving mirror trajectory defined by Eq. J4); (b) Moving mirror velocity (vertical axis) as function of the boundary position. 



III. DIAGONAL STATES: VACUUM AND THERMAL STATES 



As examples of initial field states such that the density matrix is diagonal in the Fock basis, let us consider the 
vacuum and the thermal state. It can be shown that the expected value of the energy density operator T = (Tbo(t, x)) 
can be split in T = 7~ vac + T non _ V acj where T vac is the contribution to the energy density due the vacuum part and 
^"non-vac is the non-vacuum contribution due to the real particles in the initial state of the field. Hereafter we consider 
the averages (...) taken over initial field states annihilated by B. Let us start considering the vacuum as the initial 
field state (T" non _ vac = 0). The vacuum contribution to the energy density inside the oscillating cavity can be written 
as T vac = -/(«) - /(«), where: 
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For the static situation the function R is given by R(z) — z/Lq and its first derivative is a constant R (z) = 1/Lq- 
From this equation we get the known static Casimir energy densities: T va( P = T va( ^ = — 7r/(24L§) [l|, Q and 
T™ = T^D = 7r/(48L§) [13, where the superscripts DD, NN, DN and ND means the types of boundary conditions 
considered in the calculations. 

For oscillatory laws of motion, the energy density and the force acting on the moving boundary have been considered 
in the literature. Here, we investigate the following particular trajectory x = L(t), analogous to the one proposed by 
Walker and Davies [18j |: 



t = B-L + A e-^ B - L ^ B -1 
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where we assume L(0) = L = 1, A and B are constants, with A > B so that \L\ < 1. This is a smooth and 
asymptotically static trajectory at t — > oo (see Fig. [l](a)). The mirror velocity can be relativistic (values greather 
than 0.3 of the light velocity, see Fig. [TJb)) for the parameters A = 2 and 5 = 1, considered in the presente paper. 
With these features, the energy density and the radiation force cannot be investigated via approximate methods, 
which require oscillatory motions with small amplitudes. To visualize the behavior of the energy density and the 
radiation force, we will use next the exact approach proposed by Cole-Schieve @, [HI]. 

In Fig. HJa) we plot, for both DD and NN cases and the vacuum as the intial field state, the time evolution of 
the actual force acting on the moving boundary (solid line) for each position L, whereas the dotted line shows the 
value of the static Casimir force — 7r/[24L 2 ] which would act on the boundary if it was static at the position x = L. 
In analogous manner, in Fig. [3fb) we plot, for the DN and ND cases, the time evolution of the force acting on the 
moving boundary (solid line), whereas the dotted line shows the repulsive static Casimir force 7r/[48i(t) 2 ]. We see 
that for DD and NN, and also for DN and ND, the radiation force acting on the moving boundaries are the same. In 
both cases, as expected for this law of motion, the figures show the dynamical force approaching to the static Casimir 
one in the asymptotic limit t — > oo. In Fig. [3] we see the energy density for several boundary conditions. 
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FIG. 2: (a) The quantum force acting on the moving boundary (solid line) and the attractive static Casimir force (dotted line) for the 
DD or NN cases; (b) The quantum force acting on the moving boundary (solid line) for the DN or ND cases, and the repulsive static 
Casimir force (dotted line). In both cases the vacuum was considered as the initial field state. 




(aj (b) 

FIG. 3: (a) The energy density for DD or NN cases (solid line) and also for DN or ND cases (dotted line), at the point x = Lq/2, as 
function of time; (b) The energy density for for DD or NN cases (solid line) and also for DN or ND cases (dotted line), as function of the 
normalized position x/L(t) in the cavity, for a time t = 30Lq. In both cases the vacuum was considered as the initial field state. 



Now, let us consider the thermal state as the initial field state. For this case we have (ajj<z n ') = 8 nn in{n,0) and 
(a n a n >) = (ajjajj,) = 0, where n(n,/3) = [exp(K,(n + ft)) — 1] 1 and k — 1/T. In Fig. [4] we compare the behavior 
of T n on-vac for T = 1 and T vac . We also see that T^_ VSuC =T^_ VSuC , Tj^_ VM =Tj^_ vac , and that the difference 
between TV^. vac and T^?^._ v „ c is a scale factor. 



IV. NON-DIAGONAL STATES: COHERENT AND SCHRODINGER CAT STATES 

The coherent state can be defined as an eigenstate of the annihilation operator: a n \a) = aSn n „ \ a). where a = \a\e 
and no is related to the frequency of the excited mode pj|. The Schrodinger cat state [20, l2l|, is defined as a 
superposition of coherent states | \E r ) = N (|a) + e 1 ^ \— a)), where |— a) has the same amplitude than \a) but with 
a phase shift of n, and the normalization N is given by N = (2 + 2e~ 2 ' a l 2 cosf/))^ 1 / 2 . When \a\ is as small as 2, 
| (a | — a) | 2 w [23] ■ Then we expect that the cat state behaves like two coherent states. In Fig. [5] we obtain the 
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FIG. 4: (a) The quantum forces Tnon-vac(i, L(t)) (solid line) and -30xT v ac(i, L(t)) (dotted line) for DD or NN cases, as function of 
time; (b) The quantum forces Tnon-vac (t, L(t)) (solid line) and 30 X Tvac (t, Lit)) (dotted line) for DN or ND cases, as function of time. 
In both cases the thermal bath with temperature T = 1 was considered as the initial field state. 



behavior of the force acting on the moving boundary. For \a\ = 1, we can observe differences between the non- vacuum 
forces for coherent and Schrodinger cat states. 

V. FINAL COMMENTS 

In the present paper we investigated the problem of the energy density inside a non-static cavity. Since the 
law of motion considered here presents relativistic velocities and also large displacement, we could not have used the 
usual approximate approaches found in the literature. However, the problem can be studied via an exact numerical 
approach, taking as basis the one purposed by Cole and Schieve. In the context of this method of calculation, we 
observed that the energy density and the force on the moving mirror are not affected if we change NN by DD boundary 
conditions, for vacuum or thermal initial field states. Similar behavior is observed by changing DN by ND conditions 
in the context of these diagonal states. On the other hand, the same invariance is not observed for the non-diagonal 
states investigated, for instance, the coherent and Schrodinger cat states. For these states each one of the cases (DD, 
NN, ND and DN) has a particular behavior inside the non- static cavity. This work was supported by CNPq and 
CAPES-Brazil. 
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FIG. 5: The non-vacuum part of the force (Tnon-vac(i, L(t))) acting on the moving boundary (vertical axis), as function of time. The 
solid line is the coherent case, whereas the dotted line is the Schrodinger cat case. The law of motion considered is given by Eq. Q. The 
parameters considered are: |a| = 1, 8 = <fi = 0. (a) The DD case; (b) The NN case; (c) The DN case; (d) The ND case. 
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